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Abstract. We investigate the spectral shape, the anisotropy 
of the wave vector distributions and the anisotropy of the am- 
plitudes of the magnetic fluctuations in the Earth's magne- 
tosheath within a broad range of frequencies [10^^, 10] Hz 
which corresponds to spatial scales from ^ 10 to 10^ km. 
We present the first observations of a Kolmogorov-like iner- 
tial range of Alfvenic fluctuations SB^ ^ f^^^^ in the mag- 
netosheath flanks, below the ion cyclotron frequency fd- In 
the vicinity of fd, a spectral break is observed, like in so- 
lar wind turbulence. Above the break, the energy of com- 
pressive and Alfvenic fluctuations generally follow a power 
law with a spectral index between —3 and —2. Concerning 
the anisotropy of the wave vector distribution, we observe a 
clear change in its nature in the vicinity of ion characteristic 
scales: if at MHD scales there is no evidence for a domi- 
nance of a slab (fc|| ^ k±) or 2D (fcj^ ^ fc||) turbulence, 
above the spectral break, (/ > fd, kc/ojpi > 1) the 2D tur- 
bulence dominates. This 2D turbulence is observed in six 
selected one-hour intervals among which the average proton 
/3 varies from 0.8 to 9. It is observed for both the trans- 
verse and compressive magnetic fluctuations, independently 
on the presence of linearly unstable modes at low frequencies 
or Alfven vortices at the spectral break. We then analyse the 
anisotropy of the magnetic fluctuations in a time dependent 
reference frame based on the field B and the flow velocity V 
directions. Within the range of the 2D turbulence, at scales 
[1, 30]kc/ujpi, and for any f3 we find that the magnetic fluc- 
tuations at a given frequency in the plane perpendicular to 
B have more energy along the B x V direction. This non- 
gyrotropy of the fluctuations at a fixed frequency is consis- 
tent with gyrotropic fluctuations at a given wave vector, with 
k± ^ /c||, which suffer a different Doppler shift along and 
perpendicular to V in the plane perpendicular to B. 
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1 Introduction 

In the space plasma turbulence, the presence of a mean mag- 
netic field B gives rise to anisotropics with respect to the 
field direction (|| means parallel, and _L means perpendicu- 
lar to B). There are anisotropics both in the intensities SB^ 
of the magnetic fluctuations {5B^ ^ ^-^iP 
tribution of their wave vectors k (fcj^ fcy), i.e., the energy 
distribution of the turbulent fluctuations is anisotropic in k- 
space. 

To study the anisotropy of turbulent fluctuations in space 
plasma, we chose here the Earth's magnetosheath as a lab- 
oratory. Downstream of the bow shock, the solar wind 
plasma slows down, and the plasma density, temperature 
and magnetic field increase in comparison with the so- 
lar wind plasma. The magnetosheath boundaries (bow 
shock and magnetopause) introduce an important temper- 
ature anisotropy T± > Ty, and therefore linearly unsta- 
ble waves, such as Alfven Ion Cyclotron (AIC) and mirror 
modes, are present (see the reviews by Schwartz et al., 1996; 
Lucek et al., 2005; Alexandrova, 2008). In the vicinity of the 
bow-shock, an /^^ power law spectrum is observed at fre- 
quencies below the ion cyclotron frequency, / < fd, (Cza 
jkowska et al.,,2001j ). The power law spectra ^ f^^^^, typ- 
ical of the solar wind inertial range at / < fd, have not 
been observed in the magnetosheath. However, as in the so- 
lar wind, the energy of the magnetic fluctuations f ollows a 
power law close to ~ /"'^ at frequencies / > fd (Rezeau 
[etaLj[T999l[Czaykowska et al.|[25nT] ). 

The question of the anisotropy of wave vectors in the mag- 
netosheath has been mostly addressed for dominant frequen- 
cies in the turbulent spectrum (spectral peaks), below fd, 
where linearly unstable modes are expected (Sahraoui et al. 
2004; Alexandrova et al., 2004; Schafer et al., 2005; Narita et 
al., 2006; Narita and Glassmeier, 2006; Constantinescu et al., 
2007). Instead, we are interested in permanent fluctuations in 
the magnetosheath (and not in spectral peaks) which cover a 
very broad range of frequencies (more than 5 decades), from 
frequencies well below fd to frequencies much higher than 
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These permanent fluctuations within the frequency range 
[0.35, 12.5] Hz, above fd, and for one decade of scale 
lengths around Cluster separations (~ 100 km), have been 
studied by Sahraoui et al. (2006) using the fc-filtering tech- 
nique. For a relatively short time interval in the inner mag- 
netosheath (close to the magnetopause) and for a proton beta 
(3p ^ 4, the authors show that the wave-vectors of the fluc- 
tuations are mostly perpendicular to the mean magnetic field 
B, fc_L » fcy , and that their frequency ujq in the plasma frame 
is zero. In the plane perpendicular to B, the k-distribution is 
non-gyrotropic, more intense and with a well-defined power 
law k~^/^ in a direction along the flow velocity V which was 
perpendicular to both B and the normal to the magnetopause 
for this particular case. The presence of hnearly unstable 
large scale mirror mode during the considered time interval 
makes the authors conclude that the small scale fluctuations 
with the observed dispersion properties k±_ ^ fcy andwo = 
result from a non-Unear cascade of mirror modes. 

At higher frequencies, ^ [10, 10'^] Hz, between about the 
lower hybrid frequency fih and 10 times the electron cy- 
clotron frequencie fee, the permanent fluctuations observed 
in the magnetosheath, during four intervals of several hours, 
have been studied by Mangeney et al. (2006) and Lacombe et 
al. (2006). The corresponding spatial scales, ~ [0.1, 10] km 
~ [0.3, 30] fcc/wpe (c/ujpe being the electron inertial length), 
are much smaller than the Cluster separations, and so only 
the one- spacecraft technique could be used to analyze the 
anisotropy of wave vector distributions. 

Magnetic fluctuations with k\\ » fc_L, usually called slab 
turbulence, have rapid variations of the correlation function 
along the field and weak dependence upon the perpendicu- 
lar coordinates. For the fluctuations with k±_ ^ fcy, called 
2D turbulence, the correlation function varies rapidly in the 
perpendicular plane, and there is no dependence along the 
field direction. So, measurements along different directions 
with respect to the mean field can give the information on 
the wave vector anisotropy. Under the assumption of con- 
vected turbulent fluctuations through the spacecraft (i.e., the 
phase velocity of the fluctuations is small with respect to 
the flow velocity), these measurements are possible with one 
spacecraft thanks to the variation of the mean magnetic field 
B direction with respect to the bulk flow V. While V || B, the 
spacecraft resolve fluctuations with k||B, when V _L B, the 
fluctuations with k _L B are measured. 

This idea was already used in the solar wind for study- 
ing the wave vector anisotropies of the Alfvenic fluctuations 
in the inertial range (Matthaeus et al., 1990; Bieber et al., 
1996; Saur & Bieber, 1999). The authors suppose that the 
observed turbulence is a linear superposition of two uncorre- 
cted components, slab and 2D, and both components have a 
power law energy distribution with the same spectral index 
s, 5Bl{k\\) - Aikj" and SBl{ki_) - A2fcX^ where Ai 
and A2 are the amplitudes of slab and 2D turbulent compo- 
nents, respectively. Bieber et al. (1996) propose two inde- 
pendent observational tests for distinguishing the slab com- 
ponent from the 2D component. 



The first test is based on the anisotropy of the power spec- 
tral density (PSD) of the magnetic fluctuations in the plane 
perpendicular to B, i.e., on the non-gyrotropy of the PSD at a 
given frequency in the spacecraft frame: in the case of a slab 
turbulence with k||B, all the wave vectors suffer the same 
Doppler shift depending on the angle between k and V, and 
if the spectral power is gyrotropic in the plasma frame, it 
will remain gyrotropic in the spacecraft frame; in the case 
of a 2D turbulence, with k _L B, if the PSD is gyrotropic in 
the plasma frame, it will be non-gyrotropic in the spacecraft 
frame because the Dopier shift will be different if k is per- 
pendicular to V and if k has a component along V. 

The second test reveals the dependence of the PSD at a 
fixed frequency on the angle between the plasma flow and 
the mean field Qbv (defined between and 90 degrees): For 
a PSD decreasing with k (Uke a power-law, for example), 
in the case of the slab turbulence, the PSD for a given fre- 
quency will be more intense for Q^y = 0°, and therefore, 
the PSD decreases while 6 bv increases; for the 2D turbu- 
lence the PSD will be more intense for Qbv = 90°, and so 
it increases with 6 bv- Using these tests, Bieber et al. (1996) 
have shown that the inertial range of the slow solar wind is 
dominated by a 2D turbulence; however, a small percentage 
of a slab component is present. 

Mangeney et al. (2006) proposed a model of anisotropic 
wave vector distribution without any assumption on the in- 
dependence of the two turbulence components. In their gy- 
rotropic model, the wave vector can be obUque with respect 
to the 11 and ± directions. The authors introduce a cone aper- 
ture of the angle 0kB between k and B, as a free parameter 
of the model. They assume a power law distribution of the 
total energy of the fluctuations ^ fc~*, with .s independent 
on 6kB- For a given k, the turbulent spectrum is modeled by 
one of the two typical angular distribution cos{6kB)^ for 
k nearly parallel to B and ^ sm{6kB)'^ for k nearly perpen- 
dicular to B. For these two distributions, the cone aperture 
of OkB is about 20° for /x = 10 and 7° for n = 100. The 
angle B^b can be easily represented through O bv ™d so the 
model can be tested with one- spacecraft measurements. 

An advantage of the magnetosheath with respect to the so- 
lar wind in ecliptical plane is that the angle Qbv covers the 
range from 0° to 90° within rather short time periods (one 
hour, or so) while other plasma conditions remain roughly 
the same. A comparison of the model described above with 
the observations of the total PSD of the magnetic fluctu- 
ations within the magnetosheath flanks, at frequencies be- 
tween //;, and 10/ce, shows that these fluctuations have a 
strongly anisotropic distribution of k, with O^b = (90 ± 7)° 
(Mangeney et al., 2006). Actually, this model (as well as the 
tests of Bieber et al., 1996) is valid not only for a power law 
energy distribution in k, but for any monotone dependence 
where the energy decreases with increasing k. 

Mangeney et al. (2006) have also shown that the variations 
of dB^ with Qbv for a given frequency was not consistent 
with the presence of waves with a non-negligible phase ve- 
locity Vcf,. In other words, if the observed turbulent fluctua- 
tions are a superposition of waves, their has to be much 
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smaller than the flow velocity for any wave number k. This 
is consistent with the assumption that the wave frequency l^-'o 
is vanishing: the fluctuations are due to magnetic structures 
frozen in the plasma frame. These results have been obtained 
in the magnetosheath flanks for / > 10 Hz, at electron spa- 
tial scales ~ [0.3,30]fcc/a;pe- 

In this paper we extend the study of Mangeney et al. 
(2006) to frequencies below 10 Hz, for the same time periods 
in the magnetosheath flanks. As a result, we will cover the 
largest possible scale range, from electron (~ 1 km) to MHD 
scales 10"'' km). At variance with the previous study, 
we analyse the spectral shapes and anisotropics for parallel 
(~ compressive) JSy and for transverse (~ Alfvenic) 5B± 
fluctuations independently. For Alfvenic fluctuations SB± 
we perform the first test of Bieber et al. (1996), i.e., we an- 
alyze the gyrotropy of the PSD of the magnetic fluctuations 
in the plane perpendicular to B, at a given frequency, as a 
function of 8^^. 

2 Data and methods of analysis 

For our study we use high resolution (22 vectors per second) 
magnetic field waveforms measured by the FGM instrument 
(Balogh et al., 2001). Four seconds averages of the PSD of 
the magnetic fluctuations at 27 logarithmically spaced fre- 
quencies, between 8 Hz and 4 kHz, are measured by the 
STAFF Spectrum Analyser (SA) (Cornilleau-Wehrlin et al., 
1997). Plasma parameters with a time resolution of 4 sec- 
onds are determined from HIA/CIS measurements (Reme et 
al., 2001). 

2.1 Magnetic spectra and decomposition in 6B± and 

High resolution FGM measurements allow to resolve Im- 
bulent spectra up to 10 Hz. Similar to Alexandrova et 
al. (2006), we calculate the spectra of the magnetic fluctu- 
ations in the GSE directions X, Y and Z, using the Morlet 
wavelet transform. The total power spectral density (PSD) 
is SB^if) = J2j=x.Y.z ^B^if)- The PSD of the compres- 
sive fluctuations (/) is approximated by the PSD of the 
modulus of the magnetic field. This is a good approximation 
when 5B^ ^ Bq, where Bq is the modulus of the magnetic 
field at the largest scale of the analysed data set. The PSD of 
the transverse fluctuations is therefore 

6Bl{,f) = SB\f)-SBl{f). (1) 

This approach, based on wavelet decomposition, allows the 
separation of dB± and 6B\\ with respect to a local mean field, 
i.e. to the field averaged on a neighbouring scale larger than 
the scale of the fluctuations. The lower frequency limit of 
this approach is a scale where the ordering SB"^ <^ Bq is no 
longer satisfied. 

The STAFF-SA instrument measures the spectral matrix 
{SBi{f)SBj{f)) at higher frequencies. Because of a recently 
detected error about the axes directions in the spin plane (O. 
SantoHk, 2008, private communication) we cannot separate 



parallel and perpendicular spectra at the STAFF-SA frequen- 
cies; however we present here the total PSD, the trace of the 
spectral matrix. 

2.2 Anisotropy of the k distribution 

The motion of the plasma with respect to a probe allows a ID 
analysis of the wave vector distribution along the direction of 
V, as was discussed in section 1 . The 3D wave vector power 
spectrum /(k) = I{k, OkBi^k) depends on the wave number 
k, on the angle 9kB between k and B, and on the azimuth 
of k in the plane perpendicular to B. If is the frequency 
of a wave in the plasma rest frame (wq and u are assumed to 
be positive), the Doppler shifted frequency / = a;/27r in the 
spacecraft frame is given by 

w=|wo+k-V|. (2) 
The trace of the power spectral density at this frequency is 

5B'^{u:) =Aj I{\<i)5{uj - |a;o k • V|)rfk (3) 

i.e. the sum of the contributions with different k. A is a 
normaUsation factor and 5 the Dirac function. 

The angle OkB can be considered as depending on the an- 
gle 9kv between k and V, the angle appearing in the Doppler 
shift frequency, and on the angle Qbv between B and V (see 
equation (2) of Mangeney et al., 2006). Thus, the variations 
of SB"^ with Qbv for a given lo will give information about 
/(k). As was discussed in section 1, 5B^ increases with in- 
creasing ^Bv when the fluctuations have k±_ ^ k\\ (2D tur- 
bulence) and it decreases for a slab turbulence with k\\ ^ k± 
(see Figure 6 of Mangeney et al., 2006). As a consequence, 
in the case of 2D turbulence, the spectrum of the fluctuations 
will be higher for large angles Obv than for smaU ones, and 
vice- versa for the slab geometry. 

2.3 ^B-anisotropy in the BV-frame 

To study the distribution of the PSD dBj_{f) in the plane 
perpendicular to B, i.e. the gyrotropy of the magnetic fluc- 
tuations, taking into account the direction of the flow ve- 
locity V, we shall consider the following reference frame 
(b, bv, bbv): b is the direction of the B field, bv the direc- 
tion of B X V and bbv the direction of B x (B x V). 

The definition of this frame depends on the considered 
scale (frequency). A local reference frame (defined on a 
neighbouring scale larger than the scale of the fluctuations) 
can be defined only for frequencies below the spacecraft spin 
frequency fspin = 0.25 Hz which limits the plasma mo- 
ments time resolution to 4 s. That is why, for any frequency 
/ > fspin we shall use the frame (b,bv,bbv) redefined 
every 4 s. 

In this frame, we only consider the frequencies below 
10 Hz, i.e., file FGM data (flie STAFF-SA data cannot be 
used because of the error that has to be corrected in the 
whole data set). We project the wavelet transform of Bx, 
By and Bz on the b,bv,bbv directions and we calculate 
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Fig. 1. FGM and STAFF-SA/Cluster data on 19/12/2001, 02:00- 
04:00 UT. Average spectra of the magnetic fluctuations, calculated 
using the Morlet wavelet transform of the FGM data (/ < 10 Hz). 
Solid line: for the transverse fluctuations SB±. Dashed line: for the 
compressive fluctuations . Dotted line: the total power spectral 
density, STAFF-SA data (/ > 8 Hz). The diamonds give the scales 
kc/cupi ~ krgi ~ 0.01 to 100. The vertical dotted line gives the 
average proton cyclotron frequency The shapes of the power 
laws f~^'^ and f~^'^ are shown. 



the squares of these projections 5B'^{f,t), SB^^{f,t) and 
^-^bbv (/' which are the diagonal terms of the spectral ma- 
trix in this new frame. 



3 k-distribution of 5Bj^ and SB\^ 

3.1 A case study with /3p ~ 1 

We consider an interval on the day 19/12/2001, from 02:00 to 
04:00 UT. For this interval, the mean plasma parameters are 
the following: the magnetic field B = (18±2) nT, the proton 
plasma density Np — (7 ± 1) cm""^, the proton temperature 
Tp = (120 ± 15) eV, the proton plasma beta f3p = (1.1 ± 
0.4), the ion inertial length c/u!pi — (90 ± 5) km and the ion 
Larmor radius r^j = (65 ± 15) km. The average upstream 
bow shock angle 9bn calculated with the ACE data is about 
70° (Lacombe et al., 2006). 

Figure [T] displays the average spectra of the FGM data for 
the transverse fluctuations (solid lines) and for the compres- 
sive fluctuations (dashed lines). The total PSD of the STAFF- 
SA data is the dotted line above 8 Hz. The total covered fre- 
quency range is more than six decades, from 3 • 10^"* Hz to 
300 Hz. The small vertical bars just above the abscissae-axis 
indicate the scales from A = lO"* km to 1 km corresponding 
to the Doppler shift / = V/X for 9kv = 0° and for the av- 
erage velocity V — (246 ± 25) km/s. The diamonds above 
the abscissae indicate the scales kc/uopi ~ kvgi ~ 0.01 to 
100, corresponding to the frequency / = kV/2Ti. Precisely, 
kc/ujpi = 1 appears in the spectrum at / = (0.44 ±0.05) Hz 
and krgi — 1 appears at / = (0.63 ± 0.11) Hz. 




(5R^(f = 0.52 Hz") 



20 40 60 





20 40 60 BO 



^(f = 0.2 Hz) 




20 40 60 80 



(5B^ii{f = 0.2 Hz) 




20 40 60 80 
e>„ (degrees) 



20 40 60 80 
©Bv (degrees) 



Fig. 2. FGM/Cluster data on 19/12/2001, 02:00-04:00 UT. Upper 
panels: Scatter plots of the power spectral density of the magnetic 
fluctuations at / = 3 Hz as a function of the angle between the local 
mean field and velocity, Obv- The distributions of the energy of 
Alfvenic fluctuations 5B^ is shown in the left panel, (JBjj is shown 
in the right panel. Middle and lower panels have the same format, 
but here the frequencies are respectively 0.52 Hz and 0.2 Hz. In all 
panels, the thick lines give the median value for bins 5° wide. 



We see in Figure [T] that, in the FGM frequency range, 
5B']_{f) is everywhere larger than (/), except around 

/ ~ 5 • 10^^ Hz where 5B\ ^ '^^f ' where the com- 
pressive fluctuations display a spectral break. The spectrum 
of 5B± displays a bump and a break around 0.2 Hz, that can 



be a signature of Alfven vortices (Alexandrova et al. 20061. 
Below the bump, 5B\{f) ~ /~^ **, a power law with an 
exponent close to the Kolmogorov's one —5/3 (in section 
5 we will analyse spectral shapes in more details). Above 
the bump, for kc/iOpi > 0.2, 5B\{f) and 5B'j{f) foflow a 
similar power law ^ f^"^'^ ■ This power law extends on the 
STAFF-SA frequency range up to kc/ojpi ~ 50 (kc/ujpe — 
1.3). It is quite possible that, above these scales, the dis- 
sipation of the electromagnetic turbulence starts. However, 
around / ~ 100 Hz, there is another spectral bump, which is 
due to whistler waves, identified by their right-handed polar- 
isation. The question of the turbulence dissipation is out of 
scope of the present paper and will be analysed in details in 
the future. 
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Now, we consider the anisotropy of the distribution of the 
wave vectors. Figure |2] shows the dependence of SB^ (left 
column) and (5i?^ (right column) on the angle Qbv at dif- 
ferent frequencies. The thick solid curves give the median 
values for bins 5° wide. The upper panels of Figure |2] corre- 
spond to / = 3 Hz (kc/ujpi ~ 7). The observed increase of 
5B^ and 5B^^ with Q bv can be produced only by fluctua- 
tions with k±^ ^ fcji , with phase velocities negligible with 
respect to the plasma bulk velocity, and with decreasing in- 
tensity of the fluctuations with increasing k (as was discussed 
in sections 1 and 2.2). 

At larger scales (lower frequencies) we observe the same 
tendency for / > 0.3 Hz. The middle panels of Figure |2] 
display SB^ and (5 -By as functions of Qbv for / = 0.52 Hz 
(kc/ujpi ~ 1.2): we still observe here a clear increase of SB^ 
and with Qbv- The lower panels of Figure |2j display 
SB^ and JSy at 0.2 Hz (kc/ujpi ~ 0.5), just at the spec- 
tral bump of 5B\ (see Figure jTj). (5i?y still increases with 
Q BV (in spite of a large dispersion of the data points around 
the median), while 5B\ has a flat distribution with Qbv- 
This can be due to several reasons: (i) I{k) is no longer a 
decreasing function with k, (ii) I{OkB) is more isotropic in 
the spectral bump and/or (iii) the fluctuations are not frozen 
in plasma at this scale. This spectral bump, as we have al- 
ready mentioned, can be the signature of Alfven vortices with 
k± ^ /c||, propagating slowly in the plasma frame. It can be 
also the signature of propagating AIC waves with fc|| ^ k±, 
which are unstable for the observed plasma conditions (Man- 
geney et al., 2006; Samsonov et al., 2007). However, as ex- 
plained in sections 1 and 2.2, the energy of fluctuations with 
fc|| ^ k± would decrease with increasing Qbv at a given 
frequency, while in our case 6Bj^ seems to be independent 
on Qbv- 

As we have just seen from Figure 2 (and as discussed in 
sections 1 and 2.2), the comparison of the energy level of the 
turbulent fluctuations for different Qbv at a given frequency 
gives us a good estimate of the wave vector anisotropy. We 
compare now the PSD of the fluctuations observed for large 
Q BV and for small Q bv in the whole frequency range, to es- 
timate the wave vector anisotropy in a large domain of wave 
vectors. 

The upper panel of Figure [3] displays the spectra of the 
compressive fluctuations SB'^^^^ for the 10% of the points of 
the sample with the largest Qbv Qo, = large angles, solid 
line), and <5i?|^^ for the 10% of the points with the small- 
est Qbv (sa = small angles, dashed line). At frequen- 
cies below 0.06 Hz (kc/ujpi = 0.1, indicated by a vertical 
solid line) the spectra SB'^^^^ ~ '^^f\ia- higher frequen- 
cies, / > 0.06 Hz (i.e., at smaller scales, kc/ujpi > 0.1) 
we observe ^^y^^ > SB"^^^^. This indicates that 2D tur- 
bulence dominates at such small scales. Close to 10 Hz, 
i.e. at the vicinity of the FGM Nyquist frequency, we see 
that <5-By;a — (^-Sysa — lO""' nT^/Hz, that is the sensitivity 
limit of the FGM instrument. Therefore, the observations at 
/ > 5 Hz are not physically reliable. 

The lower panel of Figure [3] displays the spectra for 
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Fig. 3. FGM/Cluster data on 19/12/2001, 02:00-04:00 UT. Aver- 
age power spectral density of the transverse magnetic fluctuations 
(upper panel) and of the compressive fluctuations (lower panel). In 
each panel, the solid line is the average spectrum for large Qbv 
angles, and the dashed line for small Qbv- The vertical dotted line 
gives the average fd, the diamonds indicate kc/upi — 0.01, 0.1 
and 1. 



the transverse fluctuations for large and small angles Qbv, 
^^±ia (solid line) and SB^^^ (dashed Une). We observe that 
^^±ia becomes larger than 6B^^^ at about the same scale of 
kc/iOpi ~ 0.1 as for compressive fluctuations. However, here 
within the spectral bump range, ^ [0.1, 0.3] Hz, we observe 
SB^i^ ~ SB^^^. This is consistent with our previous results 
that in this short frequency range the 2D turbulence model is 
not valid (cf. Figure [2]). A clear dominance of SB^^^ over 
^^±sa is '^hen observed for / > 0.3 Hz (kc/ujpi > 0.8, see a 
vertical solid line). 

These observations allow to conclude that, for Pp ~ 1, 
there is a change in the nature of the wave-vector distribu- 
tion of the magnetic fluctuations in the magnetosheath, in the 
vicinity of ion characteristic scale: if at MHD scales there is 
no clear evidence for a dominance of a slab or 2D geometry 
of the fluctuations, at ion scales (kc/ujpi > 0.1) the 2D tur- 
bulence dominates. This is valid for both the Alfvenic and 
compressive fluctuations. The large scale limit of the 2D tur- 
bulence is, however, different for Alfvenic and compressive 
fluctuations, and seems to depend on the presence of spectral 
features, as peaks or bumps. We analyse this point more in 
details by considering other cases. 
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Fig. 4. FGM/Cluster data on 17/05/2002, 08:30-09:30 UT. (a) Av- 
erage PSD of 5Bj_ (solid line), PSD of (dashed line), the ver- 
tical dotted line gives the average fd ; (b) the average spectrum of 
for large angles bv (solid line) and for small angles (dashed 
line), diamonds indicate kc/upi = 0.01, 0.1 and 1 (krgi = kc/cupi 
in this case); (c) same as (b), for SB±. 



3.2 Other case studies 

The comparison between the spectra for large Qbv and for 
small Qbv has been made during four other one-hour inter- 
vals, with different average Pp and different average shock 
angles 9bn- For the same intervals, Samsonov et al. (2007) 
display the observed proton temperature anisotropy and the 
corresponding thresholds for AIC and mirror instabilities. 

Figure |4] gives the results of the analysis for an interval 
(day 17/05/2002, 08:30-09:30 UT) for which Obn ^ 70°, 
Pp = (1.6 ± 0.3), B = (24 ± 3) nT, V = (190 ± 10) km/s, 
/ci = (0.37 ± 0.04) Hz, c/ujpi = (55 ± 2) km and rg, = 
(50 ± 5) km. 

Figure gives the average PSD of transverse (soUd line) 
and compressive fluctuations (dashed line). There is a spec- 
tral bump for the transverse fluctuations around 0.2 Hz. Be- 



low the spectral bump, 5B\{f) ^ f^^^^- For the compres- 
sive fluctuations there is a spectral bump around 0.07 Hz, 
probably made of mirror modes. Below the bump, (5i3|(/) is 
close to /^^. 

In the two other panels of Figure |4j we display the spectra 
for large and small Qbv for compressive and for transverse 
fluctuations, respectively. In Figure|4j3, at frequencies above 
the bump of 5i?y (/ > 0.1 Hz, kcjiOpi > 0.3) we observe 
SBfi, > SBla- In Figured, we observe SBl^^ 6Bl^^ 
at large scales (observed at / < 0.05 Hz, i.e. kc/cHpi < 0.1), 
but at frequencies above the bump of 6B^ (kc/cUpi > 1) we 
observe SB^^^ > SBj_^^. So, the transverse and compres- 
sive fluctuations have k± 3> fc|| at scales smaller than their 
respective spectral bumps. This confirms the conclusions of 
section 3.1. 

For an interval with /3p ~ 0.8 and 6bn — 70° (day 
16/12/2001, 08:00-09:00 UT) the average spectrum SBKf) 
displays a spectral bump around 0.5 Hz. The comparison 
between the spectra for large and small Qbv (not shown) in- 
dicates that above the bump, for kc/ujpi > 1, the transverse 
fluctuations can be described by the 2D-turbulence model. 
For the compressive fluctuations, this model is valid for a 
larger range of scales, kc/ujpi > 0.3. This confirms the re- 
sults obtained for /3p ~ 1.6, shown in Figure |4]as well as the 
conclusions of section 3.1. 

In an interval with a larger value of /3p (day 17/05/2002, 
11:00-12:00 UT, Pp ~ 4.5, 0bn ^ 73°), the analysis of the 
spectra for large and small Qbv (not shown) shows that the 
2D turbulence takes place for kc/uupi > 0.3 for the trans- 
verse fluctuations, and for kc/uipi > 0.2 for the compressive 
fluctuations. So, the 2D turbulence range of scales for the 
transverse fluctuations is wider in this case. 

Figure |5] shows the results of the analysis for an interval 
(day 16/12/2001, 05:30-06:30 UT) downstream of an obhque 
bow shock (Bbn ^ 50°), when /3p = (9 ± 3), B ^ {27 ± 
6) nT, V = (370 ± 20) km/s, /e, = (0.4 ± 0.1) Hz, c/tOpi = 
(30 ± 2) km and rgi = (60 ± 20) km. The wavenumber 
kc/iOpi = 1 appears in the spectrum at / = (2.0 ± 0.2) Hz 
and the wavenumber krgi = 1 appears at / = (1.1±0.3) Hz. 

Figure 5b shows that at low frequencies (i.e., at large 
scales, kcjujpi < 0.1) the spectrum for large angle (JSpQ 
dominates slightly at every frequencies. At smaller scales, 
kc/iOpi > 0.1, this dominance is more clear. Figure 5c 
shows that SB±ia > SB±sa as far as kc/ujpi > 0.01 (wave- 
lengths smaller than lO'' km), i.e., the dB± fluctuations can 
be described by the 2D-turbulence model at all the scales 
smaller than the Earth's radius. In this case, with large value 
of plasma beta, the 2D turbulence range of scales increases 
again, but the lower limit of 2D turbulence is not related to 
any spectral features, as was observed for smaller Pp. 

We may therefore conclude that, at frequencies above the 
spectral break in the vicinity of fd, the SB± and SB\\ fluctua- 
tions in the magnetosheath have k mostly perpendicular to B, 
and this is independent on Pp and on 9bn- In terms of spatial 
scales, this is valid for kc/ujpi > 0.1 (or > 1 when spectral 
features appear in the vicinity of kc/ojpi — 1). For high val- 
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Fig. 5. Same as Figure|4] but for FGM/Cluster data on 16/12/2001, 
05:30-06:30 UT. 



ues of (3p, the range of scales of the 2D turbulence seems to 
increase: for /3p ^ 10 the fluctuations have k± ^ fcy for 
kc/uipi > 0.01. This small scale spectral anisotropy is also 
independent on the presence of transverse and/or compres- 
sive spectral features (peaks) at larger scales. Nevertheless, 
for the moderate values of beta (/3p < 3), these spectral peaks 
appear as the lower limit of 2D turbulence. 



4 Gyrotropy of the magnetic fluctuations 

In this section we analyse the anisotropy of the amplitudes 
of magnetic fluctuations in the plane perpendicular to B. For 
this purpose, we use the coordinate frame based on B and V, 
(b, bv, bbv), as explained in section 2.3. 

For the same time interval as Figure [T] Figure |6] displays 
the ratio R ~ ^B^J^^bbv^ amplitude of the fluctuations 
along B X V over the amplitude along B x (B x V), in the 
plane perpendicular to B, at four fixed frequencies, as a func- 
tion of Qbv- 
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Fig. 6. FGM/Cluster data on 19/12/2001, 02:00-04:00 UT. Scatter 
plots of the ratio R — SB^^/SB^,,^ averaged over 4 s, as a function 
of the angle Qbv at (a) 10 Hz, (b) 3 Hz, (c) 1 Hz and (d) 0.3 Hz. 
The thick lines give the median value for bins 10° wide. 



In Figure 6a (/ = 10 Hz, kc/cupi = 23) and Figure 6b 
if — 3 Hz, kc/ujpi ~ 7) the ratio R is larger than 1 for 
©By ?J 20°. The median value decreases and reaches 1 or 
less for Qbv < 20°. A similar dependence was observed 
by Bieber et al. (1996) and Saur & Bieber (1999) at MHD 
scales in the solar wind, indicating the dominance of the 2D 
turbulence. 

At larger scales (Figure 6c, / = 1 Hz, kc/uipi — 2), in 
spite of the strong dispersion of R, the median values are 
slightly larger than 1 for any Qbv- the 2D turbulence still 
dominates. At an even larger scale, the scale of the spectral 
break (Figure 6d, / = 0.3 Hz, kc/upi — 0.7), the ratio R 
is strongly dispersed. The variation of the median does not 
correspond to a slab or 2D turbulence. That is in agreement 
with the results obtained within the spectral break frequency 
range in section 3.1 (cf Figure 2, lowest panel for SB±). 

The anisotropy of the magnetic fluctuations for the 
[10~^, 10] Hz frequency range is shown in Figure [t] with 
average spectra in the three directions b (dashed line), bv 
(solid line) and bbv (dotted line). The panels (a) to (f) cor- 
respond to increasing values of /3p for the six considered in- 
tervals. For each interval, the vertical solid bar indicates the 
scale kc/uipi ~ 1 and the dotted bar shows fd- 

In Figures |7^, b and d, for kc/ujpi > 1 we observe that 
the spectra of the components along b and along bbv are 
nearly equal, 6B^ ~ ^^bbv Figures |7^ and f, SB^ is 
larger than 6B^^: the fluctuations are more compressive for 
the largest values of (]p. All the panels of Figure|7]show that 
^^bv ^ ^^bbv fcc/wpi > 1. So, within the 2D turbulence 
range the PSD is not gyrotropic at a given frequency. 

As we have mentioned in section 1, the observed non- 
gyrotropy in the spacecraft frame can be due to the Doppler 
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Fig. 7. Average spectra of the magnetic fluctuations in the 
(b, bv, bbv)-frame, b is parallel to the B field (dashed line), bv 
is parallel to B x V (solid line), bbv is parallel to B x (B x V) 
(dotted line). For each of the 6 considered one-hour intervals a 
vertical dotted bar gives fd, a vertical solid bar gives the Doppler 
shifted wavenumber kc/uipi = 1. In each panel the shapes of the 
power laws , are indicated; in the high frequency range 

we show the f^" spectral shape, with s determined in section 5, see 
Figure 8. 



shift. Indeed, we have shown in section 3 that the wave 
vectors k are mainly perpendicular to B, i.e k lies in plane 
spanned by bv and bbv. Assuming plane 2D turbulence, the 
relation kj^ • (5B = holds and thus, the wave vectors along 
bbv (i.e., along the direction of the flow in the plane perpen- 
dicular to B, we denote such wave vectors kbbv) contribute 
to the PSD of SBfjy and the wave vectors along bv (kbv) con- 
tribute to the PSD of 5B\jin,. Even if /(k) is gyrotropic, the 
fluctuations (5i3fc„ with kbbv suffer a Doppler shift stronger 
than the fluctuations 5B\ji,y with kbv For 2D turbulence, this 
Doppler shift effect is more pronounced when Qbv reaches 
90°. 

This implies that, for a gyrotropic energy distribution, 
I{kbv) — I{kbbv), in the plane perpendicular to B, and if 
the energy decreases with fc, for example as a power law 
I{k) , the observed frequency spectrum 6B'^^{f) will 

be more intense than 6B^i^^{f). In other words, at the same 
frequency / in the spacecraft frame, we observe the fluctu- 



ations with I kbv I > I kbbv I- As the larger wave numbers 
correspond to a weaker intensity (for a monotone energy de- 
crease with fc), SB?i will be smaller than 5B?, 



bv ' 



Therefore, the non-gyrotropy of 5B^{f), observed here, 
could be due to the non-gyrotropy of the Doppler shift, and 
could be compatible with a gyrotropic distribution of /(k). 
This is confirmed by the upper panels of Figure|6] as far as k 
is mainly perpendicular to B, the Doppler shift is small and 
gyrotropic for small Qbv and we observe i? ~ 1, i.e. the 
PSD is gyrotropic; but for large Qbv, i? > 1. 

On the other hand, the ratio R{f) > 1, observed in Fig- 
ures 6a and 6b for Qbv — 90°, is also compatible with 
the non-gyrotropic k-distribution observed by Sahraoui et al. 
(2006) near the magnetopause, for Qbv — 90°. In this case 
study, the authors show that the turbulent cascade develops 
along V, perpendicular to B and n, where n is the normal 
to the magnetopause. In this geometry, the direction V is 
close to bbv. Therefore, a non-gyrotropic /(k) distribution 
with dB'^,^{k)/5B^i^^{k) > 1 is expected in the k-domain. 
This non-gyrotropy of wave vectors is then reinforced by the 
Doppler shift, and would give SB^^{f) /SB^^^{f) > 1 in the 
/-domain. 



5 Spectral shapes 

We have mentionned in section 3 that, at frequencies / < fd, 
below the spectral break, the spectra of the transverse fluc- 
tuations 5B\{f) follow a power law close to f~^^^ (see 
Figure 1 and 4a). For the six intervals of Figure 7, Fig- 
ure 8 displays compensated plots of the transverse spectra 
p/^5B\{f) (solid lines). On the low frequency side of these 
plots, we see that the compensated spectra oscillate around a 
horizontal line, in a frequency range which varies slightly 
from day to day: a power law /^^^'^ is thus a good approx- 
imation for the observations in this frequency range. In Fig- 
ures 8a to 8e, for which (3p is between 0.8 and 4.5, the Kol- 
mogoroff f^^^^ power law is observed below 0.06 or 0. 1 Hz, 
corresponding to scales kc/ujpi < 0.1 (see the vertical solid 
bar). 

In Figure 8f, for/3p ~ (9±3), the f^^^^ power law is only 
observed below 0.01 Hz, i.e., below the 2D turbulence range 
of 5B± (see Figure 5c). Above this frequency, as we see in 
Figure 7f, the spectra of all the components are close to an 
/^^ power law and the three spectra have nearly the same 
intensity, 5BlU) ~ 6Bl{f) ^ 5Bl,.^{f). This isotropy of 
the amplitudes of the turbulent fluctuations is natural to ob- 
serve in a high beta plasma, where the mean magnetic field 
does not play any important role. Within this frequency range 
the spectrum can be also formed by a superposition of AIC 
and mirror waves. For such high (3p, the mirror modes are 
more unstable than the AIC waves, but they have an impor- 
tant Alfvenic component (Genot et al., 2001), that can also 
contribute to make fluctuations more isotropic. 

Above fci, the spectra (5i?y(/) and 5B\{J) follow sim- 
ilar power laws, see Figure 7. The compensated spectra 
f'^6B^{f) with s between 2 and 3 are presented in Fig- 
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Fig. 8. Compensated spectra f^^^SBj^(f) (solid lines) and 
f"SB'^{f) with s indicated in each panel (dashed lines) for the 
same time periods as Figure 7. 



ure 8 by dashed lines. We see that the composed spectra 
p-^SB^^f) oscillate around a horizontal line in a few cases 
(Figures 7b, 7d and 7e), for different values of (3p. In Figure 
8a, the power law is steeper, s = 3. Actually, in this case 
the spectral bump is the most clearly pronounced of the six 
analyzed intervals. This bump is a signature of the Alfven 
vortices, which have their own spectrum ^ k^"^ or fc^^, de- 
pending on the vortex topology (Alexandrova, 2008). The 
superposition of the background turbulence with the coherent 
structures, like magnetic vortices, can produce the observed 
steep spectrum. In Figure 8c we do not observe any clear 
evidence for a power law spectrum in this frequency range. 



6 Summary and Discussion 

In the present paper, we have analysed six one-hour intervals 
in the middle of the terrestrial magnetosheath (at more than 
one hour from the crossing of the bow shock or of the mag- 
netopause). Precisely, we considered intervals in the magne- 
tosheath flanks: the local times for the three considered days 
are respectively 8, 17 and 18 h (Lacombe et al., 2006). The 
proton beta varies from one interval to another, (3p G [0.8, 9], 



that allows us to study the plasma turbulence in a very large 
range of plasma conditions. 

6.1 Spectral shape 

The spectral shape of the magnetic fluctuations in the mag- 
netosheath has been studied by several authors (see Alexan- 
drova, 2008). Rezeau et al. (1999) find a power law /"^i 
above fd in an interval close to the magnetopause. For inter- 
vals close to the bow shock, Czaykowska et al. (2001) find 
power laws around f^^ below fd, and /^^'^ above fd- But 
in these studies, intervals of 4 minutes have been analyzed, 
so the minimal resolved frequency is about 10^^ Hz. In the 
present paper, the length of the intervals allows to reach fre- 
quencies smaller than 10"'^ Hz. 

Here, we present, for the first time, the observations of 
a Kolmogorov-like inertial range for Alfvenic fluctuations 
5B\{f) ^ f^^^^ in the frequency range / < fd- It is 
clearly observed in five of the six studied intervals, those for 
which /3p < 5 and when Alfvenic fluctuations were dom- 
inant. Such a Kolmogorov power law is observed in the 
Alfvenic inertial range of the solar wind turbulence, below 
the spectral break in the vicinity of fd- The presence of such 
power law in the magnetosheath flanks is consistent with the 
estimations made by Alexandrova (2008): in the flanks, the 
transit time of the plasma is longer than in the subsolar re- 
gions, and it is much longer than the time of nonlinear inter- 
actions; therefore, the turbulence has enough time to become 
developed. 

In the high frequency range, / > fd, we generally observe 
5B\{f) and SB'^^{f) following similar power law f^^ with a 
spectral index s between 2 and 3, in agreement with previous 
studies. 

6.2 Wave- vector anisotropy 

We analysed here the anisotropy of wave-vector distribution 
of the magnetic fluctuations from 10^^ to 10 Hz. This fre- 
quency range corresponds to the spatial scales going from 
^ 10 to 10^ km (from electron to MHD scales). For this 
analysis we used a statistical method, based on the depen- 
dence of the observed magnetic energy at a given frequency 



on the Doppler shift for different wave vectors ( Bieber et al. 



[T996 ^, 'Horbur y et al.|[2005||Mangeney et al. ||2006| ). 

Within the inertial range of the magnetosheath turbulence 

(/ < fci, kc/ojpi < 1, kvgi < 1), we do not observe a 
clear evidence of wave-vector anisotropy. It can be related 
to the fact that linearly unstable modes, such as AIC modes 
with k mainly parallel to B and mirror modes with k mainly 
perpendicular to B, together with Alfven vortices with fcj^ ^ 
fc|| co-exist in this frequency range. 

However, above the spectral break in the vicinity of the 
ion characteristic scales (/ > fd, kc/cupi > 1, kvgi > 1 
and up to electron scales), we observe a clear evidence of 
2D turbulence with k± ^ fcy for both SB± and <5i3||, and 
independently on f3p, on the bow-shock geometry 6bn, and 
on the wave activity within the inertial range at larger scales. 
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This wave vector anisotropy seems to be a general property 
of the small scale turbulence in the Earth's magnetosheath. 

The range of wavenumbers of this 2D turbulence some- 
times goes down to fcc/wp, ~ 0.1 (or even to kc/ujpi ~ 
0.01), but usually it is limited by kc/tOpi ~ 1, while at 
kc/upi < 1 spectral features (peaks or bumps) appear. As 
we can conclude from the work of Mangeney et al. (2006), 
the largest wavenumbers of the 2D turbulence are observed 
around kc/ujpi ~ 100, where the dissipation of electromag- 
netic turbulence begins. This last conjecture must be verified 
by a deeper analysis. 

6.3 Anisotropy of magnetic fluctuations 

Analyzing the anisotropy of the amplitudes of turbulent fluc- 
tuations, we usually find that 5B\ > (5Bj| ; but for the largest 
plasma Pp, the fluctuations are more isotropic SBj^ ^ 
This is valid for both the large scale inertial range and the 
small scale 2D turbulence. The dominance of happens 
only locally in the turbulent spectrum, indicating the pres- 
ence of an unstable mirror mode. 

Concerning the gyrotropy of the amplitude of the magnetic 
fluctuations in the plane perpendicular to B, there is no uni- 
versal behavior at large scales. At smaller scales, within the 
frequency range [0.3 — 10] Hz and for any f3p, the 2D tur- 
bulence is observed to be non-gyrotropic: the energy SBf^ 
along the direction perpendicular to V and B is larger than 
the energy (5i?^j„ along the projection of V in the plane 
perpendicular to B. This non-gyrotropy might be a con- 
sequence of different Doppler shifts for fluctuations with k 
paraUel and perpendicular to V in the plane perpendicular to 
B. The non-gyrotropy at a fixed / is compatible with gy- 
rotropic fluctuations at a given k. On the other hand such a 
non-gyrotropy will be also observed if the k-distribution is 
not gyrotropic, but is aligned with the plasma flow, as was 
observed by Sahraoui et al (2006) in the vicinity of the mag- 
netopause. 
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